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ABSTRACT 


The  objective  of  this  study  is  to  examine  the  sensitivity  of 
the  low  wavenumber  portion  of  the  wall  pressure  field  under  a 
turbulent  boundary  layer  to  both  flow  and  wall  conditions.  The 
analytical  procedure  makes  use  of  a  method  for  calculating  the 
broadband  correlation  function  of  the  wall  pressure  field  in 
which  the  specific  flow  characteristics  enter  through  profiles  of 
the  mean  velocity  shear  and  turbulence  intensity  across  the 
boundary  layer.  The  basic  assumption  used  to  extend  the  results 
of  this  procedure  to  narrowband  frequency  characteristics  is 
that  there  is  a  range  of  large  spatial  separations  in  which 
Strouhal  scaling  is  not  important.  Several  specific  boundary 
layers  are  considered  including  flows  on  smooth  and  rough  walls 
and  in  zero,  favorable,  and  adverse  pressure  gradients.  Results 
are  presented  for  the  point  pressure,  the  spatial  correlation 
function,  and  the  low  wavenumber  spectral  level. 
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I.  INTRODUCTION 

When  a  turbulent  boundary  layer  develops  along  a  surface,  the  surface 
is  acted  upon  by  the  pressure  fluctuations  generated  by  the  flow.  The 
resulting  response  of  the  surface  and  the  pressure  field  it  generates  in 
many  circumstances  is  an  important  noise  source.  The  characteristics  of  the 
pressure  fluctuations  beneath  a  turbulent  boundary  layer  (TBL)  have  therefore 
been  a  subject  of  interest  for  many  years. 

Perhaps  the  simplest  conceptual  model  of  the  turbulent  flow  in  boundary 
layer  is  that  of  frozen  turbulence  in  which  a  pattern  of  random  pressure 
field  is  convected  by  the  flow  past  the  boundary  surface.  Any  pressure 
fluctuation  of  temporal  frequency  f  that  is  present  in  this  field  would  then 
have  a  spatial  wavelength  of  X^_  =  Uc/f,  where  Uc  is  the  convection  speed  of 
the  pattern.  The  associated  convective  wavenumber  is 


In  underwater  applications  where  flow  speeds  are  relatively  low  and  frequencies 

of  interest  can  be  relatively  high,  the  convective  wavelength  is  small  (and, 

conversely  the  wavenumber  is  large).  For  example  when  U  =20  kt,  X  s  0.40  in. 

c  c 

at  1000  Hz.  Because  of  their  small  wavelength,  the  noise  associated  with 
pressure  disturbances  having  spatial  scales  comparable  to  Xc  is  relatively 
easy  to  control  by  separating  the  flow  from  surfaces  of  interest  by  thin 
intermediate  layers. 

Although  measurements  of  the  spatial  scales  of  the  pressure  field  beneath 
turbulent  boundary  layers  confirm  that  the  spatial  scale  of  most  of  the 
fluctuations  is  near  the  convective  wavelength,  there  is  present  to  some  degree 
a  broad  distribution  of  wavelengths.  In  particular,  a  small  portion  of  the 
pressure  field  exists  in  disturbances  with  long  wavelengths  (i.e.,  low  wave- 
numbers)  .  When  the  scale  of  the  pressure  field  exceeds  the  convective  wave¬ 
length  by  the  Mach  number  (i.e.,  M  =  uc/c0»  where  c^  is  the  sound  speed  in  the 
fluid) ,  the  boundary  layer  noise  cannot  be  distinguished  from  an  acoustic 
pressure  disturbance  of  the  same  frequency.  The  level  of  the  pressure 
fluctuations  of  the  TBL  having  low  wavenumber  spatial  scales  therefore  sets 
a  noise  floor  for  acoustic  measurements  made  in  proximity  to  the  flow. 
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Because  the  low  wavenumber  content  of  the  wall  pressure  field  is  small 

compared  with  the  portion  near  the  convective  wavenumber,  it  has  been  difficult 

to  establish  the  level  of  the  low  wavenumber  spectrum.  Measurements  making  use 

of  wave-vector  filters  to  reject  the  high  wavenumber  components  have  been 

successful  in  setting  levels  in  small  ranges  of  wavenumber  at  moderately  low 

wavenumbers  (e.g.,  3  in  1  (Ref.  1)).  These  measurements  however  have  been 

insufficient  for  determining  the  shape  of  the  low  wavenumber  spectrum.  The 

conclusion  drawn  from  theoretical  arguments  based  on  a  planar  incompressible 

boundary  layer  of  infinite  extent  is  that  the  wavenumber  spectrum  should  behave 
2 

as  k  in  the  low  wavenumber  range.  An  analytical  model  of  the  entire  wavenumber 
spectrum  of  the  wall  pressure  fluctuations  that  embodies  this  dependence  is 
defined  in  Ref.  2. 

The  objective  of  the  present  study  is  to  examine  the  sensitivity  of  the 
low  wavenumber  portion  of  the  wall  pressure  field  under  a  TBL  to  flow  and  wall 
conditions.  The  procedure  makes  use  of  an  existing  analytical  method  for 

3 

calculating  the  broadband  spatial  correlation  function.  This  method,  which 
assumes  a  planar  flow,  is  attractive  because  it  is  formulated  explicitly  in 
terms  of  the  profiles  of  mean  velocity  and  turbulence  intensity  across  the 
boundary  layer.  The  wavenumber  spectrum  is  obtained  from  the  Fourier 
transformation  of  the  correlation  function.  Although  the  procedure  is  exact 
for  broadband  properties,  it  is  argued  that  in  the  range  of  low  wavenumber  and 
frequency  where  Strouhal  scaling  is  not  expected,  the  broadband  and  narrowband 
spatial  characteristics  should  be  nearly  the  same.  With  this  approximation 
the  procedure  is  used  to  examine  flow  with  both  zero  (i.e.,  nearly  planar)  and 
non-zero  (i.e.,  substantially  non-planar)  pressure  gradients. 

A  preliminary  study  using  the  proposed  method  is  reported  in  Ref.  4  in 
which  the  shear  contribution  of  the  wall  pressure  field  is  examined  for  the 
flow  over  a  smooth  wall  in  zero  pressure  gradient.  In  the  present  study  results 
for  both  shear  and  turbulence  contributions  are  calculated  for  four  different 
boundary  layer  flows.  Two  of  these  are  non-planar,  and  although  this  violates 
an  assumption  of  the  analysis,  they  are  included  to  evaluate  the  effect  of  a 
reasonably  large  variation  in  mean  boundary  layer  properties  on  the  results. 


-2- 
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II .  ANALYTICAL  MODEL 


A.  Definitions 


Both  temporal  and  spatial  characteristics  of  the  wall  pressure  field 
under  a  TBL  are  determined  from  the  second-order  statistics  of  the  field. 
Specifically  the  space-time  correlation  function  of  the  pressure  field  is 
defined  as  follows: 

C  (x,C,t,T)  =  E[p(x,t)p(x+5,t+T) J  (1) 

where  £  is  the  spatial  separation  vector,  t  is  the  time  delay,  and  E  is  the 
expected  value.  Under  the  assumptions  of  spatial  and  temporal  stationarity, 
the  correlation  function  becomes  only  a  function  of  separation  and  is  symmetric 
about  zero.  Furthermore  as  commonly  implemented  the  expected  value  is  taken  to 
be  the  time  average.  When  separation  time  t  is  zero,  we  obtain  the  broadband 
correlation  function 

CPP(^'0)  =  <p2>R(^}  (2) 

where  the  first  quantity  is  the  mean-square  pressure  and  R(£)  is  the  spatial 
correlation  function  normalized  to  unity  at  zero  spatial  separation. 

Narrowband  frequency  information  is  obtained  by  Fourier  transforming  the 
temporal  correlation  function.  For  example  the  cross-spectral  density  is 
defined  as  the  temporal  transform  of  the  space-time  correlation  function. 


C  (5 
PP  s 


;w)  ■  J  cpp(?'T)  e  1WT  d-r 


A  normalized  cross- spectral  density  can  be  defined  analogous  to  Eq.  2  such 
that 

C  (i,u)  =  P(u))R(5,w)  (4) 

PP 

where  R(0,oj)  =  1. 

The  spatial  Fourier  transform  of  either  Eqs.  2  or  4  yields  the  wavenumber 
spectrum  of  the  pressure  field.  For  example  when  Eq.  4  is  transformed,  the 
wavenumber- frequency  spectrum  is  obtained, 


4>  (k,ui) 


C  (Ciw)  eik’C  d? 
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where  the  integral  is  over  all  separations  on  the  surface.  Owing  to  the 
symmetry  of  the  cross-spectral  density,  the  Fourier  transform  can  be  replaced 
by  a  cosine  transform.  This  particular  function  is  useful  for  describing  the 
characteristics  of  the  wall  pressure  field. 

B.  Scaling 

Dimensionless  parameters  that  characterize  the  wall  pressure  field 

of  the  TBL  are  obtained  by  selecting  appropriate  length  and  time  scales.  For 

this  study  we  choose  for  these  scales  the  outer  parameters  of  displacement 

thickness  6^  and  mean  flow  speed  Uq.  From  these  parameters  two  length  scales 

can  be  defined,  one  being  6^  and  the  second  being  Uq/o).  The  reciprocal  of 

the  latter  is  a  wavenumber  proportional  to  the  convective  wavenumber  (k  = 

c 

w/Uc,  where  Uc  s  .6  Uq) .  For  underwater  applications  these  two  length  scales 
are  quite  disparate.  The  cross- spectral  density  from  Eq.  4  can  then  be  non- 
dimensionalized  as  follows: 

C  (5,w)U 

C  =  —EE— - -  =  P(ufi  /u  )R((4/U  ,1/6 ,)  (6) 

pp  2,  1  o  o  1 

g  61 

2 

where  q  =  l/2p  U  . 

o  o 

C.  Low  Wavenumber  Assumptions  and  Spectrum 

The  boundary  layer  is  assumed  to  be  an  incompressible  flow  actino  on 
an  infinite  plane.  Additionally  changes  along  the  flow  direction  are  neglected 
in  comparison  with  changes  across  the  boundary  layer.  This  latter  assumption 
is  compatible  with  the  previous  assumption  of  spatial  stationarity.  Under 
these  conditions,  the  theorem  of  Kraichnan  and  Phillips5  that  there  is  no  net 
force  on  the  plane  generated  by  the  TBL  is  applicable.  A  statement  of  this 
theorem  applied  to  narrow  frequency  bands  is 

00 

J  Cpp(C^)dC  =  0  (7) 

—00 


An  equivalent  statement  using  Eq.  4  is  that  the  wavenumber- frequency  spectrum 
must  vanish  at  zero  wavenumber. 
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The  low  spatial  wavenumber  levels  of  the  pressure  field  under  a  TBL 
reflect  the  information  from  the  correlation  measurements  made  at  large 
separation  distances.  At  large  separations  however,  it  is  unlikely  that  the 
Strouhal  dependence  appearing  in  Eq.  6  remains  as  an  important  parameter.  The 
standard  conceptual  argument  against  preserving  Strouhal  scaling  in  this  limit 
is  that  it  would  imply  low-frequency  correlation  of  the  wall  pressure  field  at 
very  large  separations.  Given  the  distortion  to  boundary  layers  due  to  shear, 

g 

however  this  sustained  correlation  is  not  likely  to  exist. 

The  experimental  results  for  the  scaling  of  the  large  separation  narrow- 
band  correlation  are  somewhat  conflicting.  Lack  of  Strouhal  scaling  in  a  low 

frequency/large  separation  regime  has  been  reported  in  several  experimental 

6*~8  2 
studies.  Additionally,  the  analytical  model  developed  by  Chase'2  is  consistent 

7 

with  the  data  trends  of  Bull.  Measurements  reported  in  Refs.  9  and  10  however 
do  not  indicate  a  lack  of  Strouhal  scaling  at  low  frequency.  One  explanation 
by  Blake10  of  this  discrepancy  with  the  results  of  Bull  is  that  the  latter  data 
are  distorted  by  the  bandwidth  of  a  filter  in  the  measurements  system. 

In  the  present  study  we  assume  that  there  is  a  frequency/separation 
regime  in  which  Strouhal  scaling  fails  and  we  examine  the  implications  to  the 
(narrowband)  wavenumber-frequency  spectrum.  Under  this  assumption,  the  cross- 
spectral  density  in  Eq.  6  is  given  by 

?  >  «L 

S  (5,n;0))  =  P(U)6  /U  )R(5/6.  ;  (8) 

pp  1  O  11 

n  >  nL 

where  £  and  n  are  the  separation  distances  in  the  flow  and  transverse  directions, 
respectively,  and  the  subscript  L  indicates  the  separation  beyond  which  Strouhal 
scaling  fails.  It  is  understood  that  the  frequency  range  in  Eq.  8  is  also 
limited  by  the  assumed  lack  of  Strouhal  scaling. 

When  both  sides  of  Eq.  8  are  integrated  over  frequency  the  simplification 
introduced  by  the  assumed  independence  of  Strouhal  scaling  becomes  apparent. 

The  result  is  the  definition  of  the  broadband  spatial  correlation  function  given 
in  Eq.  2.  Consequently,  under  the  assumption  of  no  Strouhal  scaling  the  spatial 
characteristics  of  the  wall  pressure  field  at  large  separation  become  independent 
of  the  frequency  bandwidth. 
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The  wavenumber  frequency  spectrum  is  obtained  from  the  two-dimensional 
Fourier  transform  of  the  cross-spectral  density  function.  The  integrals  are 
over  all  separation  distances  and  involve  contributions  from  both  large  and 
small  separation  distances.  In  particular  the  narrowband  wavenumber  spectrum 
contains  contributions  from  the  narrowband  correlation  function  at  small 
separations  which  involves  Strouhal  scaling  (i.e.,  Eq.  6).  Invoking  Eq.  8,  we 
obtain  the  following  expression  for  the  narrowband  wavenumber  spectrum: 


^(k1,k2;w)  =  4P  (to) 


r«L  "l 

If 


R(£  ,  r),to/UQ)  cosk^£cosk2nd£dri 


0  0 


4  f 


J  j  R(C/n)cosk15cosk2nd£:dn 


CL  \ 


(9) 


The  Kraichnan-Phillips  theorem  expressed  in  Eq.  7  states  that  this  spectrum 
must  be  zero  at  k^  =  0.  When  applied  to  Eq.  9  as  well  as  to  its  broadband 
counterpart,  this  theorem  gives  the  following  identity: 


L 


n 


L 


R(£,n»to/U  )cosk  nd£dn  = 
o  z 


R(£,n)cosk2nd£dn 


(10) 


This  equality  states  that  the  average  over  £  of  the  narrowband  and  broadband 
correlation  functions  at  small  separations  are  equal.  By  Eq.  9  each  of  these 
contributions  equals  the  negative  of  the  ^-average  over  large  separations. 

The  relationship  in  Eq.  10  strictly  applies  to  zero  wavenumber  (k^) .  We 
assume  that  a  similar  relationship  between  integrals  of  the  narrowband  and 
broadband  correlation  function  exists  at  small  but  non-zero  wavenumber.  Under 
this  assumption  the  low  wavenumber  integral  can  be  replaced  by  its  broadband 
counterpart,  giving  the  following  approximate  result  for  the  low  wavenumber- 
frequency  spectrum: 


Lim  $  (k. ,k  ,-w) 


4P  (w) 


ff 


R(£,  n)  cosk  £cosk  r)d£dr) 


(11) 
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This  result  is  used  in  this  report  to  obtain  the  low  wavenumber  frequency 
spectrum  from  a  model  of  the  broadband  correlation  function  R(£,n) . 

D.  Spatial  Correlation  Function 

The  broadband  spatial  correlation  function  in  Eq.  11  is  evaluated 
using  the  analytical  procedure  developed  in  Ref.  3  by  Meecham  and  Tavis 
which  assumes  the  turbulent  boundary  layer  to  be  a  planar  incompressible  flow. 
The  formalism  and  numerical  results  to  facilitate  implementation  of  the 
procedure  are  also  given  in  Ref.  3;  consequently,  only  an  outline  of  the 
method  is  given  here. 


The  pressure  field  in  the  incompressible  boundary  layer  is  governed  by 
the  following  equation: 


V 


2 

P  = 


-P 

o 


(12) 


where  the  subscripted  variables  are  Cartesian  tensors  with  indices  1  and  3 
referring  to  the  flow  and  normal  directions,  respectively,  and  the  underlined 
variables  are  vectors.  The  source  term  on  the  right-hand  side  is  given  by 

Ha3(r')  =  2U(z')6alUg  +  uaug  -  <uaug>  .  (13) 


where  6. .  =  1  is  the  Kronecker  delta  function  (6. .  =  1,  i=j;  6. ,  =0,  i^j) , 

13  ijij 

U  and  u  are  the  mean  and  fluctuating  flow  speeds,  respectively,  and  z  is  the 
coordinate  in  the  normal  direction.  This  equation  which  holds  in  the  semi¬ 
infinite  space  above  a  rigid  wall  is  solved  in  terms  of  an  infinite  space 
Green's  function  by  replacing  the  wall  by  an  image  source  field  and  a  drag 
force.  The  latter  quantity  is  shown  in  Ref.  3  to  be  negligible.  Under  these 
conditions  the  solution  for  the  pressure  field  in  the  boundary  layer  is  given 
by 


P(r) 


3  H  (r  ) dr 
a3  -o  -o 

3r  3r 
oa  o3 


(14) 


where  integral  is  over  all  space. 

The  normalized  correlation  function  defined  by  Eqs.  1  and  2  is  obtained 
by  time-averaging  the  product  of  the  pressures  at  two  spatial  locations, 
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Cpp(r,r')  =  <p(r)p(r')> 


(15) 


When  Eq.  14  is  substituted  into  Eq.  15,  the  correlation  is  expressed  in  terms 
of  a  2-fold  volume  (i.e.,  6-fold  linear)  integral  with  multiple  point  velocity 
correlations  appearing  in  the  integrand.  Simplifications  are  introduced  by 
using  the  assumption  of  isotropic  turbulence  to  model  these  correlation 
functions.  The  scale  of  the  turbulence  is  taken  to  depend  on  the  distance 
from  the  wall.  Furthermore  under  the  assumption  of  a  Gaussian  distribution 
function,  fourth-order  velocity  correlations  are  expressed  in  terms  of  second- 
order  correlations.  With  this  assumption  and  extensive  manipulations  four  of 
the  integrals  can  be  evaluated  analytically  in  terms  of  relatively  simple 
functions  of  the  velocity  correlations.  The  remaining  two-fold  integrals  have 
to  be  evaluated  numerically. 


Two  contributions  to  the  correlation  of  the  surface  pressure  are  obtained. 
The  first  which  depends  on  the  correlation  of  the  fluctuating  velocity  components 
is  given  by 


(2)  Cdd)(T1) 

R  2  (n)  =  -EE- — 


w  uu 

-Jj  J  ac3[u,(|c3|)/u0]4/H  J 


Q(n,A,C3)<3A 


(16) 


-n 


where  Q  is  an  explicit  function  given  in  Ref.  3  in  terms  of  several  tabulated 
functions,  M  is  the  scaling  function,  and  5 ,  A  are  transformed  separation 
coordinates.  This  contribution  has  been  shown  to  be  small  in  general.  Results 
for  the  contribution  of  this  term  to  the  mean  square  pressure  are  presented  in 
Section  IV.  The  second  contribution  to  the  correlation  depends  on  the  product 
of  the  mean  velocity  gradient  and  the  turbulence  velocity.  This  contribution 
can  be  expressed  as  follows: 


rpp,({-"’ 


0  L  J  0 


C3)pdp 


(17) 


where  N  is  defined  in  terms  of  several  tabulated  functions,  and  p,  are 
transformed  separation  coordinates. 

As  indicated  in  Eqs.  16  and  17,  the  spatial  correlation  of  the  wall 
pressure  is  given  in  terms  of  two-fold  integrals  that  have  to  be  evaluated 
numerically.  The  inner  integral  is  a  "universal"  function  in  that  it  is 
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independent  of  specific  boundary  layer  parameters.  In  this  formulation 
variations  between  boundary  layers  introduced  by  roughness  or  pressure 
gradients  are  obtained  through  the  profiles  of  mean  velocity  and  of 
turbulence  intensity  that  appear  in  the  outer  integral. 


-9- 
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III.  NUMERICAL  EVALUATION 

A.  Universal  Functions 

The  functions  Q  and  N  in  the  inner  integrand  of  Eqs.  16  and  17  are 
defined  in  Ref.  3  in  terms  of  explicit  integrals  of  functions  of  the  velocity 
correlations.  A  smooth  representation  of  the  velocity  correlation  function 
which  is  required  to  evaluate  these  integrals  is  used  in  Ref.  3.  This 
reference  also  contains  all  the  functions  required  to  calculate  Q  and  N  both 
in  graphical  and  in  tabular  form.  In  the  present  study  the  tabulated  data  are 
used  to  calculate  the  required  integrals,  these  data  being  interpolated  using 
third-order  Lagrange  and  logarithmic  (for  large  arguments  of  monotonically 
decreasing  functions)  interpolation  schemes.  The  inner  integral  in  Eqs.  16 
and  17  is  evaluated  using  Simpson's  rule. 

B.  Turbulent  Boundary  Layer  Data 
1.  Data  Sources 

The  mean  flow  data  required  in  the  analytical  calculation  are 

the  profiles  of  turbulence  intensity  and  mean  velocity  shear  across  the  boundary 

layer.  For  the  purpose  of  the  present  study  these  data  are  required  for  flows 

over  smooth  and  rough  flat  walls  with  pressure  gradients  that  range  from 

favorable  to  adverse.  Sets  of  measurements  that  include  both  turbulence  and 

mean  velocity  data  for  these  conditions  are  limited  to  a  small  number  of  studies. 

Data  for  zero  pressure  gradient  flows  on  a  smooth  wall  from  a  number  of  studies 

are  reported  in  Hinze's  book.11  Of  particular  note  are  the  measurements  by 

12 

Klebanoff  of  the  turbulence  intensity  close  to  the  wall  (i.e.,  z/<5  <  .02). 
Measurements  on  a  smooth  wall  of  both  mean  velocity  and  the  longitudinal 
component  of  turbulence  intensity  for  zero,  adverse  and  favorable  pressure 
gradients  are  reported  by  Schloemer.13  Data  for  a  flow  with  zero  pressure 
gradient  on  smooth  and  rough  walls  are  reported  by  Blake1^  and  by  Burton.1^ 
Finally,  measurements  on  both  smooth  and  rough  walls  in  a  range  of  pressure 
gradients  is  given  by  Burton.15  This  latter  reference  is  particularly  useful 
in  that  the  profiles  of  turbulence  intensity  in  several  directions  are  given 
for  these  boundary  layers. 

Specific  parameters  of  the  boundary  layer  flows  selected  for  this  study 
are  given  in  Table  I.  The  corresponding  turbulence  intensity  and  mean 
velocity  data  are  discussed  in  subsequent  sections. 
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2.  Data  Curve  Fit 


Turbulence  Intensity 


The  analytical  model  which  involves  a  numerical  integration 
across  the  boundary  layer  requires  the  profile  of  turbulence  intensity.  This 
profile  is  defined  by  the  experimental  meas cements  at  a  relatively  small 


number  of  distances  from  the  boundary.  Furthermore  in  those  data  sets  that 


include  more  than  one  component  of  the  turbulence  intensity,  the  measurement 


locations  can  differ  for  each  component.  When  more  than  one  turbulence 
velocity  component  is  reported,  we  use  the  following  definition  for  the 
turbulence  intensity: 


•'  ■  * 


where  ui  is  the  root-mean- square  turbulence  velocity  in  the  i  direction  and 


N  components  are  measured  (i.e.,  N  =  1,  2,  or  3). 


In  order  to  use  the  discrete  measurements  in  the  analytical  model,  we 


must  fit  a  continuous  curve  through  the  data  points.  Once  such  a  curve  has 
been  defined,  the  data  can  be  interpolated  at  any  location  between  measurement 
points.  Spline  interpolation  is  used  for  this  purpose.  A  cubic  spline  is  a 
third-order  polynominal  that  fits  through  two  data  points.  The  set  of  splines 
selected  to  span  the  data  is  then  constrained  to  result  in  a  function  having 
minimum  curvature.  This  constraint  results  in  adjacent  polynomials  being 
joined  continuously  with  continuous  first  and  second  derivatives  (see  Ref.  16 
for  further  discussion  and  FORTRAN  program) . 


The  generation  of  an  interpolating  function  for  the  turbulence  intensity 
data  occurs  in  two  steps.  First  separate  spline  fits  are  made  to  the  data  for 
each  component  of  the  turbulence  intensity,  this  being  necessary  for  those 
cases  where  the  measurement  locations  differ  for  each  component.  The  specific 
sets  of  data  used  are  shown  in  Fig.  1.  The  number  of  data  points  for  each 
component  typically  varies  from  10  to  15,  and  the  data  range  between  z/6  =  .05 
and  z/S  s  1.0.  The  second  step  is  to  use  these  functions  to  calculate  the 
turbulence  intensity  according  to  Eq.  1.  Except  for  the  data  of  Ref.  12  (see 
Fig.  2)  there  are  no  measurement  available  to  define  the  variation  of  intensity 
near  the  wall.  It  is  therefore  assumed  that  the  behavior  in  all  boundary  layers 


vv  w  .  '  .  ■ .  • ...' 
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is  similar  to  that  of  the  zero  pressure  gradient  flow  along  a  smooth  wall. 
Specifically  it  is  assumed  that  for  z/d  <  0.05  the  dependence  on  distance 
from  the  wall  is  as  follows: 


I 

u 


I 

u 


- (z/.05) 2 


(19) 


where  z  =  z/d.  The  measured  data  away  from  the  wall  along  with  the  above 

assumption  near  the  wall  are  used  to  generate  the  turbulence  intensity  profile 

-4 

at  approximately  20  points  spanning  the  range  of  z/d  from  2  x  10  to  1.2. 


b.  Mean  Velocity  Shear 

With  the  possible  exception  of  the  flow  in  an  adverse  pressure 

gradient  the  mean  velocity  profile  near  the  wall  is  well  represented  by  a 

17 

universal  logarithmic  distribution  of  the  form, 


=  2.5  Ln(z/k  )  +  B  (20) 

*  s 

u 


where  k  is  the  equivalent  sand  roughness  and  B  is  a  parameter  depending  on 
s 

the  hydraulic  roughness  of  the  wall  (see  Fig.  3) .  When  the  wall  is  hydraulically 
smooth  (i.e.,  roughness  height  less  than  the  viscous  sublayer  thickness),  B  is 
given  by 

I  u*k  \ 

B  =  2.5  Ln  \—~J  +  C  (21) 

where  v  is  the  kinematic  viscosity  of  the  fluid,  and  C  is  given  by  various 
sources  to  be  in  the  range  of  4.5  -  5.5.  Substituting  of  Eq.  21  into  Eq.  20 
results  in  the  familiar  form  of  the  universal  distribution  for  smooth  walls, 
namely, 

-  2.5  Ln  Hi  +  c  (22) 

rr  v 


For  completely  rough  walls  B  takes  on  a  constant  value  of  8.5. 

The  mean  velocity  shear  of  all  boundary  layers  having  a  velocity  profile 
given  by  Eq.  20  is  the  same,  namely, 
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It  is  obvious  that  this  relationship  does  not  apply  in  the  immediate  proximity 
of  the  wall.  In  the  case  of  a  smooth  wall,  there  is  a  constant  stress  layer 
where  the  mean  shear  can  be  related  to  the  shear  velocity  and  the  fluid 
viscosity.  Since  the  mean  shear  profile  is  multiplied  in  the  analytical  model 
by  the  profile  of  turbulence  intensity  which  vanishes  at  the  wall,  a  complete 
description  of  the  mean  shear  is  not  required. 

The  universal  velocity  distribution  in  Eq.  20  does  not  extend  to  the 
outer  reaches  of  the  boundary  layer.  To  account  for  the  outer  region,  the 
expression  in  Eq.  20  can  be  made  non-universal  by  adding  the  correction  term 
provided  by  Cole's  law  of  the  wake.  This  term  can  be  adjusted  to  match  the 
deviations  from  the  universal  law  that  are  encountered  with  rough  walls  and 
with  non-zero  pressure  gradients.  However,  in  the  outer  region  of  the  boundary 
layer  both  the  turbulence  intensity  and  the  mean  shear  (see  Eq.  23)  become 
small.  Furthermore,  since  these  two  profiles  are  multiplied  together  in  the 
calculation,  it  is  not  necessary  to  include  the  shear  contribution  from  the 
law  of  the  wake  term. 

Equation  23  has  been  used  for  the  mean  velocity  shear  for  all  the 
boundary  layers  in  zero  and  favorable  pressure  gradients.  The  flow  in  the 
adverse  pressure  gradient  reported  in  Ref.  15  however  is  not  found  to  follow  a 
universal  logarithmic  distribution.  For  this  flow  the  velocity  distribution  is 
reasonably  well  described  by  a  power  law, 

u/Uq  =  (z/<5)1/n  (24) 

with  n  s  1.50  for  the  flow  with  a  free  stream  velocity  of  U  =  100  ft/s  (see 

o 

Fig.  4) .  Using  this  form  for  the  velocity  profile  along  with  the  definition 
of  the  displacement  thickness,  we  find  that 

V5  -  ir  a  °-4  1251 

The  mean  velocity  shear  in  this  boundary  layer  is  given  by 


This  shear  has  a  weaker  dependence  on  distance  from  the  wall  than  that  of 
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3.  Data  Interface  with  Numerical  Integration 

The  numerical  calculation  of  the  shear-turbulence  contribution 

to  the  correlation  function  requires  the  product  of  turbulence  intensity  and 

mean  velocity  shear.  For  convenience  this  product  is  evaluated  at  approximately 

-4 

20  points  from  z/S  =  2  x  10  to  z/6  =1.2  using  the  spline  fit  to  the  turbulence 
data  discussed  in  Section  B.2  and  the  analytical  expression  for  the  mean  shear 
discussed  in  Section  B.2.  A  cubic  spline  fit  to  the  product  is  made  which  is 
then  called  by  the  numerical  integration  routine. 

The  results  for  the  turbulence  intensity,  mean  velocity  shear,  and  the 
product  for  flow  with  zero  pressure  gradient  along  smooth  walls  (Fig.  5);  zero 
pressure  gradient  along  rough  walls  (Fig.  6);  favorable  pressure  gradient  along 
smooth  walls  (Fig.  7) ;  and  adverse  pressure  gradient  along  smooth  walls 
(Fig.  8).  The  outer  integral  in  Eqs.  16  and  17  is  evaluated  using  Simpson's 
rule.  The  number  of  terms  required  for  convergence  of  these  integrals  is  found 
to  depend  on  the  spatial  separation. 

C.  Wavenumber  Spectrum 

The  wavenumber  spectrum  is  given  by  the  two-fold  cosine  transformation 

of  the  correlation  function  (see  Eq.  11) .  These  integrals  are  evaluated  using 
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Filon's  integration  formula.  The  shear  noise  contribution  to  the  correlation 
differs  from  the  purely  turbulence  contribution  in  that  it  depends  on  the 
direction  of  spatial  separation  (e.g.,  flow  or  transverse).  Because  of  this, 
the  two-fold  transformation  of  the  shear  noise  contribution  becomes  a  product 
of  the  one-dimensional  transforms  in  each  direction,  that  is, 


$(kl'k2) 


4P  (w) 


/ 


R(C)cosk1Cd? 


/ 


r (n) cosk^hdn 


(27) 


The  narrowband  spectrum  level  in  Eq.  11  is  assumed  to  be  given  by  the  point 
pressure  spectrum  (see  for  example  Ref.  19). 
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IV.  RESULTS 


In  this  section  results  are  presented  for  the  point  pressure,  the  spatial 
correlation  functions,  and  the  low  wavenumber  level  of  the  four  flows  discussed 
in  Section  III.  These  boundary  layers  (viz.,  zero  pressure  gradient  on  smooth 
and  rough  walls,  and  both  favorable  and  adverse  pressure  gradient  flows  on 
smooth  walls)  are  selected  to  be  representative  of  the  range  of  variation  of 
mean  and  turbulence  profiles  encountered  in  practice.  Because  of  its  highly 
non-planar  flow,  the  boundary  layer  in  the  adverse  pressure  gradient  is 
inconsistent  with  the  basic  assumptions  of  the  analysis.  It  is  of  interest 
however  to  see  whether  the  calculated  results  bear  any  resemblence  with 
experimentally  determined  characteristics. 

The  root-mean- square  (RMS)  pressure  normalized  to  the  dynamic  head  is 
given  by 

P/q  =  IC  (0)1*  (28) 

1  2 

where  q  =  -r  PU  ,  and  the  "bar"  indicates  RMS  value.  If  we  normalized  the 
2  o 

pressure  to  the  wall  shear  stress  then  the  following  relationship  is  found: 


[C  (0)J 
PP 


(29) 


Measurements  of  this  quantity  by  a  number  of  investigators  in  boundary  layers 
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with  nominally  zero  pressure  gradient  are  found  to  be  in  the  range, 

P/t  =  3  ±  1  (30) 

w 

Values  of  this  ratio  have  been  calculated  for  the  four  flows  examined  in  this 
study.  Both  the  shear  and  turbulence  contributions  to  the  point  pressure  have 
been  calculated.  A  summary  of  these  results  is  given  in  Table  II  where  the 
total  pressure  and  the  shear  contribution  are  compared  with  experimentally 
determined  values.  For  all  four  boundary  layers  the  turbulence  term  contributes 
at  most  11  percent  to  the  total  RMS  point  pressure.  This  is  consistent  with 
the  findings  of  Ref.  3.  With  the  exception  of  the  flow  under  the  adverse 
pressure  gradient,  the  calculated  results  are  within  the  range  given  by  Eq.  30 
and  are  in  reasonable  agreement  with  the  specific  experimental  determination. 
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The  normalized  spatial  correlation  function  in  both  the  flow  and 
transverse  directions  is  plotted  in  Figs.  9  thru  12  as  a  function  of  separation 
distance  normalized  to  displacement  thickness.  Based  on  the  results  of  the 
point  pressure  calculations,  only  the  shear  contribution  to  the  correlation 
functions  has  been  evaluated.  The  results  for  the  flow  over  the  smooth  wall 
in  zero  pressure  gradient  shown  in  Fig.  9  indicate  a  positive  correlation  in 
the  transverse  direction  for  all  separations  and  a  correlation  in  the  flow 
direction  that  changes  sign  at  a  normalized  separation  of  approximately  3.5. 

In  general  the  magnitude  of  the  correlation  in  the  transverse  direction  exceeds 
that  in  the  flow  direction,  this  being  consistent  with  the  broadband 
measurements  of  Ref.  9  and  21.  The  results  for  the  correlation  in  the  flow 
direction  display  a  long  region  of  negative  correlation  that  slowly  returns  to 
zero  at  large  separation.  This  result  in  general  gives  a  somewhat  negative 
average  to  the  correlation  (and  by  Parseval's  theorem,  to  the  level  at  zero 
wavenumber)  in  the  flow  direction.  This  result  violates  the  Kraichman-Phillips 
theorem  (i.e.,  Eq.  7).  Since  there  is  no  experimental  information  available  in 
the  range  of  very  low  wavenumbers  (i.e.,  in  and  somewhat  above  the  acoustic 
range)  and  since  the  effects  of  this  negative  average  are  less  important  at 
somewhat  higher  wavenumbers,  this  aspect  of  the  results  is  not  thought  to 
adversely  effect  the  comparison  among  flows.  The  most  likely  reason  for  the 
long  region  of  negative  correlation  is  the  error  introduced  by  implementing 
the  procedure  to  calculate  the  correlation  function  via  the  tabulated  values 
given  in  Ref.  3.  Specifically,  the  behavior  of  the  correlation  function  at 
large  separation  depends  on  the  values  of  several  functions  at  large  argument. 
Interpolation  of  these  functions  between  the  tabulated  values  is  therefore  a 
likely  source  of  error. 

In  Fig.  10  the  correlation  functions  are  shown  for  the  boundary  layer 
flow  in  zero  pressure  gradient  over  a  rough  wall.  These  correlation  functions 
are  minimally  different  from  those  of  the  smooth  wall  flow  shown  on  Fig.  9. 

The  principle  difference  is  found  in  a  larger  point  pressure  for  the  rough 
wall  flow  (see  Table  II).  The  normalized  correlation  functions  shown  in 
Fig.  11  for  the  flow  over  the  smooth  wall  in  the  favorable  pressure  gradient 
are  also  similar  to  those  obtained  for  the  zero  gradient  flow. 
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Figure  12  shows  the  results  for  the  correlation  functions  in  the  boundary 
layer  on  a  smooth  wall  in  an  adverse  pressure  gradient.  As  discussed  above 
calculations  of  this  non-planar  flow  are  not  expected  to  agree  well  with 
measurements.  The  correlation  in  the  flow  direction  has  a  deep  region  of 
negative  correlation  compared  with  the  results  obtained  for  the  previous  flows. 
Additionally,  even  the  transverse  correlation  is  somewhat  negative  at  large 
separations . 

Results  for  the  low  wavenumber-frequency  spectrum  are  presented  on 
Figs.  13  thru  16  for  the  four  flows.  These  results  are  calculated  via  Eq.  27 
at  a  Strouhal  number  (S  =  of  unity.  For  each  flow  the  normalized 

spectrum  is  plotted  against  dimensional  wavenumber  for  a  range  of  low  wave- 
numbers  that  is  above  the  acoustic  wavenumber  (u)/co)  for  frequencies  below 
approximately  800  Hz.  Also  presented  on  these  figures  is  the  zero  wavenumber 
level  calculated  for  these  flows  using  a  modified  Corcos  spectrum.  At  a 
Strouhal  number  of  unity  for  these  flows  this  spectrum  is  essentially 
independent  of  wavenumber  with  the  normalized  level  being  given  by 


$>(0,0,u))Uc 

2^3 

q 


?  /  0  \2/  *  \< 

i  I  _£  ]  (  u  ) 

,  \U  /  \ U  / 

3  'o'  'o' 


-3  2  -2 

S  J(i+cp 


(31) 


where  the  ratio  of  convection  to  external  velocity  is  taken  to  be  0.6  and 
the  constants  are  =  .08  and  C3  =  .55. 

On  Fig.  13  the  wavenumber  spectrum  is  shown  for  the  zero  pressure  gradient 
boundary  layer  on  a  smooth  wall.  The  increasing  slope  at  low  wavenumbers  is  a 
result  of  the  somewhat  negative  value  at  zero  wavenumber.  At  the  higher  wave- 
number  range  the  level  is  approximately  12  dB  above  that  of  the  modified 
Corcos  level.  The  level  for  the  zero  pressure  gradient  flow  over  the  rough 
wall  shown  on  Fig.  14  is  found  to  be  approximately  16  dB  above  the  smooth  wall 
result  and  approximately  20  dB  higher  than  the  Corcos  spectrum.  The  level  of 
the  spectrum  calculated  for  the  favorable  pressure  gradient  flow  shown  on 
Fig.  15  is  approximately  27  dB  higher  than  that  of  the  zero  pressure  gradient 
flow  along  the  smooth  boundary,  and  it  is  substantially  higher  than  that  derived 
from  the  modified  Corcos  spectrum.  It  should  be  noted  that  the  latter  spectrum 
has  been  defined  by  adjusting  the  Corcos  spectrum  to  agree  at  low  wavenumber 
with  available  measurements,  these  resulting  from  flows  in  nominally  zero 
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pressure  gradient  on  smooth  walls.  It  is  therefore  not  unreasonable  to 


expect  a  poor  comparison  with  the  modified  Corcos  spectrum  for  the  flows  with 
non-zero  pressure  gradient.  Figure  16  shows  the  results  for  the  flow  in  an 
adverse  pressure  gradient. 
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V.  CONCLUSIONS 

The  objective  of  this  study  has  been  to  explore  the  utility  of  an 
approximate  procedure  for  evaluating  the  effects  of  flow  non-ideality  on  the 
narrowband  characteristics  of  the  wall  pressure  field,  specifically  on  the 
wavenumber- frequency  spectrum  at  low  wavenumbers.  The  procedure  which  is 
"exact"  for  broadband  results  but  approximate  for  narrowband  spectra  has  been 
implemented  using  partial  results  available  in  tabular  form  in  Ref.  3. 

Conclusions  drawn  from  this  study  are  the  following: 

1.  The  RMS  point  pressure  normalized  to  the  wall  shear  stress  that  is 
obtained  for  three  of  the  flows  is  in  reasonable  agreement  with  the  experimentally 
determined  value.  The  value  obtained  for  the  flow  in  a  substantially  adverse 
pressure  gradient  is  considerably  larger  than  that  measured. 

2.  The  levels  of  the  low  wavenumber  spectra  are  somewhat  higher  than 
the  levels  obtained  from  the  modified  Corcos  spectrum  for  the  two  zero  pressure 
flows.  Substantial  differences  are  obtained  between  these  two  spectra  for  the 
flows  in  non-zero  pressure  gradients. 

3.  The  calculated  spectra  at  very  low  wavenumbers  are  contaminated  by 
errors  introduced  by  the  use  and  interpolation  of  the  tabulated  functions. 

This  finding  is  detrimental  to  a  more  general  usage  of  the  method. 
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TABLE  I 

BOUNDARY  LAYER  PARAMETERS 


Flow 

U  (ft/s) 
o 

6  (ft) 

*1 

(ft) 

u*/U 

o 

k  (a) 
g 

(ft) 

^1  dp 
q  dx 

Smooth  Wall/ 

Zero  Pressure 
Gradient14 

164 

.142 

2.3 

x  10~2 

.034 

0 

0 

Rough  Wall/ 

Zero  Pressure 
Gradient14 

164 

.200 

4.6 

x  10-2 

.055 

7.7  x 

io-3 

0 

Smooth  Wall/ 
Favorable  Pressure 
Gradient15 

164 

.039(b) 

4.9 

x  10"3 

.049 

0 

-9.2  x  10”4 

Smooth  Wall/ 

Adverse  Pressure 
Gradient15 

100 

.180(C) 

7.2 

x  10-2 

.022 

0 

19 

(a)  Arithmetic  mean  roughness  height 

(b)  Assumed  to  be  8  6^ 

(c)  See  Eq.  25 
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TABLE  II 

CALCULATED  RATIO  OF  RMS  SURFACE  PRESSURE  TO  WALL  SHEAR  STRESS 


P/l 

w 

Flow 

u*/U 

C  (0) 

PP 

Calculated 

Total 

c(1) 

PP 

Measured 

Smooth  Wall 

Vp  =  0 

0.034 

7.2  x  10-5 

4.1 

3.7 

3.6 

Rough  Wall 

VP  =  0 

0.055 

4.7  x  10”4 

3.8 

3.6 

2.9 

Smooth  Wall 
Favorable  VP 

0.049 

1.5  x  10“4 

2.7 

2.5 

2.1 

Smooth  Wall 
Adverse  VP 

0.022 

4.5  x  10"5 

68 

67 

8 
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(a)  Variation  of  turbulence  kinetic  energy  across  the  boundary  layer. 

14 

(b)  Variation  of  turbulence  intensity  near  the  wall. 


Mean  velocity  profiles  with  inner  parameter  scaling.  Favorable 
gradients.13 


Profiles  of  normalized  mean  velocity  shear  ( - ) ,  turbulence  intensity  ( - ) ,  and 

their  product  ( - )  across  the  boundary  layer.  (Smooth  wall;  zero  pressure  gradient.) 


Profiles  of  normalized  mean  velocity  shear  ( - ) ,  turbulence  intensity  ( - ) ,  and  their 

product  ( - )  across  the  boundary  layer.  (Smooth  wall;  adverse  pressure  gradient.) 


Normalized  spatial  correlation  function  in  the  flow  ( - )  and  transverse  ( - — )  directions 

(Smooth  wall;  favorable  pressure  gradient.) 


Normalized  spatial  correlation  function  in  the  flow  ( - )  and  transverse  ( — - )  directions 

(Smooth  wall;  adverse  pressure  gradient.) 


Fig.  16  -  Normalized  wavenumber  spectrum  at  a  Strouhal  number  of  unity  calculated  using  the  present 

procedure  ( - )  and  the  modified  Corcos  model  ( — - )  .  (Smooth  wall;  adverse  pressure 

gradient. ) 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES,  INC. 


REFERENCES 


1  T.  M.  Farabee  and  F.  E.  Geib,  Jr.,  Measurement  of  Boundary  Layer  Pressure 
Fields  with  an  Array  of  Pressure  Transducers  in  a  Subsonic  Flow, 

David  Taylor  Naval  Ship  R&D  Center  Rept.  76-0031  (March  1976) . 


2  D.  M.  Chase,  "Modeling  the  Wave  Vector-frequency  Spectrum  of  Turbulent 
Boundary  Lyer  Wall  Pressure",  J.  Sound  Vib.  70  (1980)  29-67. 

3  W.  C.  Meecham  and  M.  T.  Tavis,  Theoretical  Pressure  Correlation  Functions  in 
Turbulent  Boundary  Layers,  Rept.  SAMSO-TR-78-42  prepared  for  Air  Force 
Weapons  Laboratory  (June  1984).  Also,  Phys.  Fluids  22  (1980)  1119-1131. 

)  4  j.  E.  Cole,  III,  Spectrum  of  the  Turbulent  Boundary  Layer  Wall  Pressure 

Fluctuations  at  low  Wavenumber,  CAA  Rept.  U-861-320.1  prepared  for 
Office  of  Naval  Research  (October  1982) . 


5  R.  H.  Kraichnan,  "Pressure  Fluctuations  in  Turbulent  Flow  Over  a  Flat  Plate", 

J.  Acoust.  Soc.  Am.  2J3  (1956)  378-390.  Also,  see  discussion  in  Ref.  2. 

6  M.  J.  Fisher  and  P.  O.  A.  L.  Davies,  "Correlation  Measurements  in  a  Non-frozen 

Pattern  of  Turbulence",  J.  Fluid  Mech.  JL8  (1964)  97. 

7  M.  K.  Bull,  "Wall-pressure  Fluctuations  Associated  with  Subsonic  Turbulent 

Boundary  Lyaer  Flow",  J.  Fluid  Mech.  28  (1967)  719-754. 

8  J.  M.  Clinch,  "Measurements  of  the  Wall  Pressure  Field  at  the  Surface  of  a 

Smooth-walled  Pipe  Containing  Turbulent  Water  Flow",  J.  Sound  Vib.  9 
(1969)  398-419. 


9  W.  W.  Willmarth  and  C.  E.  Wooldridge,  "Measurements  of  the  Fluctuating  Pressure 
at  the  Wall  Beneath  a  Thick  Turbulent  Boundary  Layer",  J.  Fluid  Mech  1_4 
(1962)  187. 

10  W.  K.  Blake,  Turbulent  Boundary  Layer  Wall  Pressure  Fluctuations  on  Smooth 
and  Rough  Walls,  MIT  Acoustics  and  Vibrations  Laboratory  Rept.  70208-1 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES.  INC. 


11  J.  0.  Hinze,  Turbulence,  2nd  ed. ,  {McGraw-Hill  Book  Co.,  New  York,  1975), 

626-642. 

12  P.  S.  Klebanoff,  National  Advisory  Committee  on  Aeronautics  Tech  Note 

No.  3178  (1954). 

13  H.  H.  Schloemer,  "Effects  of  Pressure  Gradients  on  Turbulent-Boundary 

Layer  Wall-Pressure  Fluctuations",  J.  Acoust.  Soc.  Am.  42^  (1967)  93-113. 

14  T.  E.  Burton,  On  the  Generation  of  Wall  Pressure  Fluctuations  for  Turbulent 

Boundary  Layers  Over  Rough  Walls,  MIT  Acoustics  and  Vibration  Laboratory 
Rept.  70208-4  (March  1971). 

15  T.  E.  Burton,  Wall  Pressure  Fluctuations  at  Smooth  and  Rough  Surfaces  Under 

Turbulent  Boundary  Layers  with  Favorable  and  Adverse  Pressure  Gradients, 
MIT  Acoustics  and  Vibration  Laboratory  Rept.  70208-9  (June  1973) . 

16  G.  E.  Forsythe,  et  al..  Computer  Methods  for  Mathematical  Computations, 

(Prentice-Hall,  Inc.,  Englewood  Cliffs,  NJ,  1975)  71-79. 

17  H.  Schlichting,  Boundary  Layer  Theory,  6th  ed.,  (McGraw-Hill  Book  Co., 

New  York,  1968)  580-583. 

18  M.  Abramowitz  and  I.  A.  Stegun,  Handbook  of  Mathematical  Formulas, 

(NBS  Applied  Math.  Series  55,  1964)  890. 

19  J.  O.  Hinze,  op  cit.,  p.  674. 

20  W.  K.  Blake,  Aero-hydroacoustics  for  Ships,  David  Taylor  Naval  Ship  R&D 

Center  Rept.  DTNSRDC-84/010  (June  1984)  733. 

21  M.  K.  Bull,  Properties  of  the  Fluctuating  Wall-pressure  Field  of  a  Turbulent 

Boundary  Layer,  AASE  Rept.  234  (March  1963) . 


-40- 


END 


FILMED 

10-85 

DTIC 


